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ABSTRACT 

Aims. A super-solar fluorine abundance was observed in the carbon-enhanced metal-poor (CEMP) star HE 1305+0132 ([F/Fe] = 
+2.90, [Fe/H] = -2.5). We propose that this observation can be explained using a binary model that involve mass transfer from an 
asymptotic giant branch (AGB) star companion and, based on this model, we predict F abundances in CEMP stars in general. We 
discuss wether F can be used to discriminate between the formation histories of most CEMP stars: via binary mass transfer or from 
the ejecta of fast-rotating massive stars. 

Methods. We compute AGB yields using different stellar evolution and nucleosynthesis codes to evaluate stellar model uncertainties. 
We use a simple dilution model to determine the factor by which the AGB yields should be diluted to match the abundances observed 
in HE 1305+0132. We further employ a binary population synthesis tool to estimate the probability of F-rich CEMP stars. 
Results. The abundances observed in HE 1305+0132 can be explained if this star accreted 3-11% of the mass lost by its former AGB 
companion. The primary AGB star should have dredged-up at least 0.2 M Q of material from its He-rich region into the convective 
envelope via third dredge-up, which corresponds to AGB models of Z 0.0001 and mass =i 2M Q . Many AGB model uncertainties, 
such as the treatment of convective borders and mass loss, require further investigation. We find that in the binary scenario most 
CEMP stars should also be FEMP stars, that is, have [F/Fe] > +1, while fast-rotating massive stars do not appear to produce fluorine. 
We conclude that fluorine is a signature of low-mass AGB pollution in CEMP stars, together with elements associated with the slow 
neutron-capture process. 

Key words. Stars: individual: HE 1305+0132 - Stars: AGB and post-AGB - Stars: binaries - Stars: abundances - nuclear reactions, 
nucleosynthesis 



1. Introduction 

Carbon-enhanced metal-poor (CEMP) stars are chemically 
peculiar objects, which represent 10-20% o f all halo stars 
dBeers & Christliebl 120051; ICohen et al.l 120051; ILucatello etail 
120061) . Most of CEMP stars exhibit radial velocity vari- 
ations, which imply the presence of a binary companion 
( ILucatello et all l2005h . A significant fra ction of CEMP stars 
(-70-80%, according to lAoki et alj|2007l CEMP-s) also exhibit 
enhancements in heavy elements such as Ba and Pb, which are 
produced by slow neutron ca ptures (s process) in asymptotic gi- 
ant branch (AGB) stars (e.g.. lGallino et al.1 fl998). One scenario 
to explain the abundance patterns in CEMP stars is therefore 
mass transfer from a former AGB companion in which the car- 



bon and heavy n e utron-capture elements were produced (e.g., 
Ilvans et all 120051 : [Thompso n et al.l 120081) . However, a certain 
fraction of CEMP stars, which have typically [Fe/Hfl< -2.7, 
exhibit low or no neutron-capture element abundances (CEMP- 
no). These stars might have formed instead from materi al ejected 
by rapidly rotating massive stars ( Meynet et al.l 120061) or faint 
Type II supernovae dUmeda & Nomotol2005l) . At extremely low 
metallicities, [Fe/H] < -4, giant CEMP stars could have en- 
riched themselves in carbon via a "dual core flash" - where mix- 
ing of protons during the core helium flash induces a hydrogen 
flash - while in the early phases of AGB stars of masses < 1.5 
M Q and [Fe/H] < -2.3, a "dual shell flash" ma y occur, where 
protons are ingested into the convective pulse (Fujimoto et al. 



* The first two authors have contributed equally to this paper. 
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Fig. 1. C, N, and F yields from AGB models of different masses 
and [Fe/H]= -2.3. Note that the F yields are plotted on a differ- 
ent scale. 



1990 t iHollowell et all 1 19901: iFuiimoto etaflkOOOt IPTcardi et alj 
2004tlCristallo et al.ll2007l). 



Schul er et al .] d2007l) derived a super-solar fluorine abun- 
dance of A(Fj3 = +4.96 ± 0.21 for the halo star HE 1305+0132, 
which corresponds to [F/Fe] = +2.9. This is the most Fe- 
deficient star, [Fe/H] = -2.5 ± 0.5, for which the fluorine abun- 
dance has been measured to date. HE 1305+0132 also exhibits 
overabundances of C and N ([C/Fe] = +2.68 ± 0.51;[N/Fe] = 
+ 1.6 ± 0.46) and an O abundance typical of halo stars ([O/Fe] 
= +0.50 ± 0.22) . Lines of Ba and Sr are observed in its spectra 
dGoswamill2005h . which place HE 1305+0132 in the group of 
CEMP-s stars. 

Fluorine abundan ces were first determined in AGB stars by 
iJorissen et al. I (119921) . Enhancements of up to 30 times the so- 
lar value were reported, demonstrating that these stars produce 
fluorine. Observations of post-AGB stars and planetary nebulae, 
the progeny of AG B stars, confirm that these objects are a lso en- 
riched in fluorine (Werner et al. 2005; Zhan g & Liull2005l) . Type 
II supernovae (IWoosley & Haxtonlll988l) and W olf-Rayet stars 
during helium burning dMevnet & Arno uld 2000) have been the- 
oretically identified as F production sites, but they are not obser- 
vationally confirmed. In contrast to that observed in CEMP stars 
and in particular HE 1305+0132, Type II supernovae typically 
produce more O tha n C, a part from a n arrow range of initial 
mass around 80 M (Wooslev et al. 2002) or in the case of faint 
supernovae at Z = dUmeda & Nomotol 12005). while models 
of Wolf-Rayet stars show that fluorine production in these stars 
scales with stellar metallicity an d decreases when rapid rotation 
is included dPalacios et al.ll2005l) . 

The aims of this paper are to discuss the yields of C, N, and F 
from AGB stars at the metallicities relevant to CEMP stars (Sect. 
2); to determine if the abundances measured for HE 1305+0132 
can be explained using the AGB binary scenario (Sect. 3); and to 
evaluate the consequences of fluorine production in AGB stars 
on the CEMP stellar population (Sect. 4). We evaluate AGB 
modelling uncertainties related to different physics prescriptions 
and nuclear reaction rates in Sect. 5 and present our conclusions 
in Sect. 6. 



2. C, N, and F from low metallicity AGB stars 

Fluorine can be produced in AGB stars via the 
18 0(p,Q') 15 N(Q', y) 19 F reaction chain during thermal pulses 
(TPs) associated with the periodic activation of the He-burning 



shell. While lg O is provided by the 14 N(a,y) 18 F03 + v) 18 O chain, 
protons are produced via 14 N(n,/;>) I4 C reactions, and neutrons 
originate from 13 C(a, «) 16 reactions. The 13 C and 14 N nuclei 
are mixed down into the convective He shell from the ashes of 
the H-burning shell. After each TP quenches, third dredge-up 
(TDU) may occur, which trasports 19 F to the convective enve- 
lope, together with 12 C produced by partial He burning in the 
TPs. 

Yi elds for C, N and F at fFe/Hl = -2.3 from the m od- 
els of iLugaro etail (120041) and lKarakas & Lattanziol d2007l) are 
shown in Fig. Q] The reader is referred to these papers for de- 
tails of the computational methods. The profiles of the C and the 
F yields as a function of the initial stellar mass closely follow 
each other. This is because the TDU carries primary 12 C to the 
stellar envelope. This is converted into primary I3 C and 14 N in 
the H-burning ashes, whose abundances drive the synthesis of F. 
Hence, fluorine p roduction in AGB st ars of low metallicity is of 
a primary nature (Lugaro et al. 2004). When the stellar mass is 
higher than 3 M , proton captures at the base of the convec- 
tive envelope (hot bottom burning, HBB) lead to the conversion 
of C into N and to the destruction of F. Hence, we expect that 
the hig h F abundance ([F/Fe] « +4) predicted by Mevnet et al. 
(120061) in the envelope of a low-metallicity rotating 7 M model 
at the beginning of the AGB phase will be completely destroyed 
by HBB during the subsequent AGB phase. 

Since the 12 C/ 13 C rati o in HE 1305+0 132 is observed close 
to its equilibrium value dGoswam H l2005l) and N is mildly en- 
hanced, proton captures may have occurred inside this star, or in 
the AGB companion, due to non-convective-mixing processes 
that occurr below the inner boundary of the convective enve- 
lope. These can be du e, for example, to efficien t activation of 
thermohaline mixing (Char bonnel & Zahnl 120071) . An alterna- 
tive explanation may be the occurrence of the dual shell flash 
in the AGB star. The high observed F abundance indicates that 
F cannot have been significantly destroyed by proton captures 
via the 19 F(p, a) 16 reaction, thus constraining the temperature 
reached by any extra mixing process to less than 20 million 
de grees, within reactio n rate uncertainties (see, e.g, Fig. 2 and 3 
of lArnould etaill 19991) . 



3. Estimate of binary transfer parameters 

In Fig.[2]we show the abundances of F and C+N,0 relative to H, 
observed in the star and compare these values to the abundance 
ratios in the material lost by AGB stars according to the models 
presented in Fig. Q] To reproduce the observed abundances we 
employ two free parameters: (1) the initial mass of the polluting 
AGB star, which determines the chemical composition of the ac- 
creted material, and (2) the amount by which the accreted mate- 
rial is diluted into the envelope of the polluted star. In Fig.|2]the 
composition that is produced by the mixing of material before 
accretion with material from the AGB star lies on a straight line 
connecting the two components; the amount of dilution is repre- 
sented by the position of the point on the mixing line. To explain 
the measured abundances, we require an AGB initial mass of be- 
tween 1.7 and 2.3 M Q , and dilution of the accreted material by a 
factor of between six and nine. 

If we assume that the mass of HE 1305+0132 is 0.8 M , 
which is the typical mass of halo stars, and calculate the evo- 
lutionary track of a star of this mass and metallicity Z = 10~ 4 
(using the stellar evolution code STARS, see, e.g., iPols et all 



A(Element) = log N(Element)+12. 



3 We prefer to use C+N rather than C because our models do not take 
into account the possibility that some C is converted into N. 
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1 1995b . we find that the observed e ffective temperature T - 
4.46 + 0.10 kK dSchuler et alj|2007l) indicates that the star is a 
giant with a convective envelope. When the envelope reaches 
its maximum depth, the outermost 60% of the mass of the star 
is convectiveQ With these assumptions, we find that the star 
should have accreted 0.05-0. 12 M from a former AGB compan- 
ion. Given that the total mass lost during AGB evolution is in the 
range 1.0-1.5 M , this corresponds to the accretion of 3-1 1% of 
the mass lost by the AGB star. 

We implicitly assume that the composition of the accreted 
material can be represented well by the average composition of 
the total material ejected by the AGB star. In reality, however, 
the binary orbit is altered during mass transfer, so the compo- 
sition of the accreted material varies because the surface AGB 
composition varies. Moreover, the evolution of the AGB star it- 
self may be altered by the presence of a binary companion. Since 
there are no models presently available that describe both binary 
and detailed AGB evolution simultaneously, we have defaulted 
to using the AGB average composition. 

In principle, it is possible to constrain the initial orbital pe- 
riod of the binary system. The system must have been suffi- 
ciently wide to allow the primary star to evolve into an AGB 
star before filling its Roche lobe; this requirement sets a lower 
limit to the initial orbital period. On the other hand, the efficiency 
of wind accretion decreases with the distance between the two 
stars, which implies that the required amount of accreted mate- 
rial may set an upper limit on the in itial orbital period. Using the 
binary evolution code described in llzzard et"al] d2006l) . we find 
that the range of initial orbital periods that corresponds to an ac- 
cretion efficiency > 3% in a system with initial masses 0.8 and 
2 M is between 7 and 27 years. 



4. Fluorine abundances in CEMP stars from binary 
systems 

With the binary evolution code described in llzzard et al. (2006), 
we sim ulated a population of binaries assuming the Kroupa et al. 
(1993) initial mass function, a uniform distribu tion of initial 
mass ratios, and the iDuq uennov & Ma yor] dl99ll) initial period 
distribution. 

We find that, of all turn-off and giant CEMP stars formed 
by mass transfer from an AGB companion, 8 1 % are expected to 
show an enhanced fluorine abundance of [F/Fe] > +1, 12% have 
[F/Fe] > +2, and 0.005% have [F/Fe] >+3. This can be under- 
stood qualitatively from the results presented in Fig. Q] because 
AGB stars that produce carbon also produce fluorine. To date, 
fluorine enhancements have not been reported for other CEMP 
stars because high-resolution spectra in the 2.3-2.4 //m band are 
required, which are not yet available. We predict that fluorine 
should be found in most CEMP stars that formed via the AGB 
binary scenario, thus representing a tracer of low-mass AGB pol- 
lution in addition to s-process element enhancements. 

The fluorine abundance in the particular case of 
HE 1305+0132 appears to be exceptionally high, since we 
expect only 0.04% of turn-off and giant CEMP stars to have 
A( 19 F)> +4.75, when considering the observational error bar. 
However, we emphasize that many uncertainties play a role in 
this estimate, both in the assumed distribution functions (which 
are reasonable for stars in the solar neighbourhood, but not 



4 The mixed fraction might be as m uch as 90% if thermohaline mix- 
ing is effective ( StanclifFe et al. 2007), even though this may be debated 
in particular reg arding the possible counter-effect of g ravitational set- 
tling (see, e.g.. lThompson e t al. 2008; Aoki et al. 2008). 
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Fig. 2. Abundances by number of F and C+N with respect to 
H as observed in H E 1305+0132 with 1 <x errors, taken from 
ISchuler et al.l {2007, hatched ellipsoid), and computed in the av- 
erage material lost via the winds by AGB stars at [Fe/H] = -2.3, 
with labels indicating the initial masses. The initial composition 
is scaled from solar. The gray area is the region of the plot where 
the F and C+N abundances from the AGB companion should lie 
to ensure that the observed abundances are reproduced after di- 
lution. The darker gray area represents the region covered by our 
models with masses in the range indicated by the asterisks. 



necessarily for halo stars), and in t he assumed wind accr etion 
efficiency. The latter is based on the lBondi & Hovkj l 19441) pre- 



script i on, whereas hydr odynamical simulations ( Theuns et al.l 
1996; Nagae et al. 2004) predict typically lower accretion 
efficiencies. These uncertainties deserve more attention in a 
follow-up study (Izzard et al., in preparation). Moreover, the 
theoretical understanding of F production in AGB stars is itself 
still affected by many uncertainties, as discussed below. 



5. Uncertainties in AGB models 

To evaluate AGB model uncertainties we discuss a set of mod- 
els of 2 M and [Fe/H] = -2.3, computed using different physics 
and nuclear reaction rate assumptions (Table [T). The first four 
models in Table Q] are com puted using the codes described in 
iKarakas & Lattanziol(l2007l) . Model 1 is that used in the previous 
section for comparison with HE 1305+0132. Model 2 includes a 
region in which protons from the envelope are mixed down into 



Table 1. Number of TPs with TDU, total TDU mass, and C, N, 
and F yields (all in M ) for different AGB models of 2 M at 
[Fe/H]= -2.3. 
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Karakas" 
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3.7E-02 


7.2E-05 
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1.7E-05 
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Cristallo r 
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iLugaro et alJ <2004h a nd Karak as & Lattanziol i2007h . iKarakas etaD 
( 2008). ICristallo et alJ d2007h . ICristallcl ( 120060 . 
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the top layer of the He- and C-rich intershell (the region between 
the H and He shells) at the end of each TDU. Proton captures on 
12 C generate a "pocket" rich in 13 C, th e main neutro n source for 
the i-process in these stars, and 14 N dHerw ig 2005). Extra 15 N 
is produced in the 13 C- 14 N pocket, which is then converted to 
19 F in the following convective TP (see also lGorielv & Mowlavil 
2000). The introduction of a 13 C- 14 N pocket of 0.002 M Q in- 
creases the 19 F yield by 60%, while the C and N yields are unaf- 
fected. In Model 3 of Table[T]we consider a model computed us- 
ing the upper limit (UL) of the 18 F(a, p) 21 Ne rate. This increases 
the 19 F yield by a f actor of ^2.5 (see details and discussion in 
iKarakas et aT]|2008l) . When both a 13 C- 14 N pocket and the up- 
per limit of the 18 F(a, p) 2l Ne reaction rate is used (Model 4 of 
Table Q]), the 19 F yield increases by a factor of three. The other 
main nuclear uncertainties originate from the 14 C (a, y) 18 Q and 
the 19 F(g, p) 22 Ne reaction rates (see discus sion in Lugaro et alj 
|2004|) . For the latter, the latest evaluation bv ltfgalde et al J (12008) 
needs to be tested in AGB models. 

All models presented do not take account of the effects in- 
duced by carbon enhancement on the opacities of the cool exter- 
nal layers of AGB stars. When C/O > 1, C-bearing molecules, 
most notably C2 and CN, increase the opacity of the external 
layers, causing th e envelope to expand and the star to become 
larger and cooler dMarigoll2002l) . Models of AGB stars of low 
mass and metallicity with C- and N-enhanced low-temperature 
opacities have been calculated using the Fr ascati Raphson 
Newton evolutionary cod e (FRANEC) code dStraniero et alj 
2006; Cristallo et al. 2007). In these models the mass-loss rate 
strongly increases with respect to models in which opacities are 
always calculated using the initial Z = 10~ 4 solar-scaled compo- 
sition. The resulting yields (Model 5 of Table [1} are w 5 times 
smaller than in the Karakas models. For comparison, the results 
obtained with the same code, but usin g opacities ca lculated for 
the initial solar-scaled composition (Cristallol l2006l) are also re- 
ported in Table Q] (Model 6), and are in good agreement with 
the Karakas model, in spite of the different choic es of mass-loss 
rate a nd treatment of the convective borders (see lStraniero et ail 
120061 for details). 

It is evident that further work is required to address the un- 
certainties in the AGB fluorine yields at low metallicities. In par- 
ticular, the inclusion of low-temperature carbon-enhanced opac- 
ities in the Karakas models (Karakas, Wood & Campbell, in 
preparation) will provide an independent comp arison to the re - 
sults obtained by the FRANEC code and by iMarigol d2002l) . 
Since a clear dependence of the mass-loss rate on the metallic- 
ity has still not be identified, different mass-loss prescriptions 
should be tested (Cristallo et al., in preparation). Finally, we 
note that the possible occurrence of the dual shell flash at the 
beginning of the AGB phase may also affec t fluorine production 
and n eeds to be investigated in detail (e.g, Camp bell & Lattanziol 
2008). 



6. Discussion and conclusions 

We have shown that it is possible to reproduce the C and F abun- 
dances observed in the CEMP-s star HE 1305+0132 via binary 
mass transfer from a companion by accretion of 3-11% of the 
mass lost by the primary star during its AGB phase. The AGB 
star should have dredged-up at least ^0.2 M of its intershell ma- 
terial into the convective envelope by means of the TDU. While 
rapidly rotating massive stars produce enough carbon and nitro- 
gen t o form CEMP stars , they do not appear to produce fluo- 
rine dPalacios et alj|200"5l) . The binary formation scenario is thus 
favoured in the case of HE 1305+0132. In general, we predict 



that most CEMP stars formed by mass transfer from an AGB 
companion should also be FEMP stars, i.e., have [F/Fe] > +1. 
Hence, fluorine appears to be a useful discriminant between the 
different scenarios proposed for the formation of CEMP stars. 

During the preparation of this manuscript, another halo ob- 
ject highly enri ched in fluorine wa s discovered, the planetary 
nebula BoBn 1 dOtsuka et al.ll2008l) . The metallicity, as well as 
all the C and N abundances observed in this object are the same, 
within errors, as those of HE 1305+0132, which suggests that 
BoBn 1 has a close connection to CEMP stars, perhaps repre- 
senting the evolutionary outcome of a CEMP star. On the other 
hand, the derived F abundance in this object is roughly a factor of 
three higher than that obtained for HE 1305+0132. This observa- 
tion can be explained via the binary scenario only if we consider 
the F yields we computed including the 13 C- 14 N pocket or the 
upper limit of the 8 F(ar, p) 21 Ne reaction. Following this indica- 
tion, we multiplied the F yields by a factor of three in our stellar 
population model and calculated a probability of 12% for CEMP 
stars to have A( 19 F)> +4.75. This result provides us with a possi- 
bility to alleviate the problem of the extremely small probability 
of the high F abundance assessed for HE 1305+0132. 
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